Laser-Induced Entangled Giant Quasi-Molecules in Optical Lattices 
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We analyze atom-atom interactions in optical lattices due to a laser-induced long-range inter- 
atomic force which prevails over the usual London-van der-Waals forces. This force, which can be 
generated by an intense laser field at a wavelength longer than that of the lattice-generating laser, 
is shown to bind pairs of cold atoms trapped at different lattice sites, and cause their translational 
quantum correlations (spatial entanglement). 

- - . PACS numbers: 34.10.+x,33.80.-b,42.50.Ct42.50.Fx 

The best-known laser-induced force acting upon or between atoms is the " dipole force" , resulting from the action 

■ of the gradient of the laser field, which has been extensively used to confine cold atoms in optical traps Q . This force 
has also been exploited in optical lattices, created by intersecting laser beams, which can trap cold atoms in periodic 
potential wells with typical depth of neV and inter- well separations on the order of the laser wavelength [|| . Recently, 
the mechanical effect of another laser-induced force has come to the fore, that of dipole-dipole interactions generated 

' by off-resonant laser beams with small gradient (nearly plane waves) ||. In the near-zone region of multiple, mutually- 
, detuned and appropriately oriented laser beams, such dipole-dipole interactions yield a gravity-like interatomic 1/r 
attraction Q|. In contrast to the previously considered near-resonant dipole-dipole interaction || via real photon 
exchange between a ground-state atom and its excited partner, wherein their cooperative decay plays a role, the far- 
off resonant dipole-dipole interaction arises from virtual (nearly decay-free) photon exchange between two ground-state 
atoms, which acquire small induced dipole-moments j3|. 

Here we put forward a new effect of the off-resonant, laser-induced dipole-dipole interaction: its ability to bind 
', together pairs of ground-state atoms localized at adjacent sites of an optical lattice. The resulting novel "giant" 
(submicron-size) quasimolecules will remain bound even if the lattice potential is switched off. An important feature 
of such quasimolecules are their translational quantum correlations or entanglement, i.e., nonfactorizability of the 
two-atom spatial wavefunction. The predicted effects should be contrasted with entanglement between internal states 
of closely separated atoms residing at the same site of a densely occupied optical lattice, via the essentially different 
near-resonant dipole-dipole interaction J5|. 

The retarded dipole-dipole interaction potential between atoms A and B, induced by an off-resonant circularly- 

■ polarized, plane- wave laser is 

Vab — — — — Fe(kr); Fg(kr) = cos(fcr cos 0) 
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{cos(fcr) + fcr sin(fcr)}(l — 3 cos 2 0) t (1 + cos 2 9) cos(kr) ^ 
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Here 9 is the angle between the interatomic axis and the wave- vector k — lu/c of the laser, I is the laser intensity and 
a = ^jgr^g) is the atomic dynamic polarizability, with d and loa being the atomic dipole moment and transition 
frequency, respectively, assuming that both atoms (A, B) have the same polarizability. In a 3D optical lattice, the Vab 
potential in Eq. (1) can be attractive for a particular relative position of two tightly-bound atoms whilst repulsive for 
another. Yet the lattice allows us to find a direction of propagation of the additional laser beam such that Vab (Ax, 9) 
is attractive for two atoms located at adjacent lattice sites, ki,r ~ 7r, along certain lattice axes, thus making the 
quasibinding effect observable, as illustrated below. 

For pair binding at adjacent lattice sites to have a distinct spectral signature, the binding potential depth |Vab(t' — 
7t//cl)| should substantially exceed two parameters: (a) the linewidth for the off-resonant absorption of the additional 
laser, nTS, where T is the resonant transition linewidth and the saturation parameter S — 4 gl^. r a ; with Q being 
the laser Rabi frequency and 5 = lu — lua its detuning, should satisfy S « 1, thus keeping the interacting atoms 
unexcited, as required by the validity of Eq. (1) (b) the heating energy TjTheat by the additional laser given by 

H Thoat ~ EjiT^2y laS " > /\ Vab\, where Er — h 2 k 2 /2m is the recoil energy of an atom with mass m, and the inelastic 
Rayleigh-scattering rate T^"y las ' ) satisfies Q ?irp"° las ' ) /|VAB| ~ fhD(k ■ r)/\Fg(kr)\, f^uik ■ f) being the Lamb-Dicke 
factor of the system, which reduces to (kr) 2 <C 1 in the near zone, and Fg(kr) is given by (1). Conditions (a) and (b) 
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can be satisfied in the near-zone regime of Vab, at r ~ 7r/fci <C 7r/fc, such that Fg(kr) 1 and /ld(A; • r) <C 1. This 
requires the additional laser to have a much longer wavelength than the lattice period, k <C kj,. 

We proceed with the general theoretical framework that allows one to treat atom-atom scattering in an optical lattice 
for any interatomic potential FAs(r), which is assumed to be a function of the interatomic separation r = | — rg| 
only. The Hamiltonian of the two interacting atoms is H = Hq + Vab , where 

h 2 h 2 

Ua(b) being the periodic potential of the lattice acting on atom A(B). 

The formation of a bound or quasi-bound pair is revealed by the analysis of the transition matrix 

T = V A B(l-GVA B y 1 (3) 

G being the unperturbed Green function (associated with Ho). Let D be the determinant of the matrix (1 — GVab)- 
When Re(Z?) = 0, a bound state (resonance with zero width) may be formed if its energy Eb exceeds the total band 
energy of the two unperturbed atoms. A resonance may occur if Re(_D) = corresponds to an energy E r within the 
unperturbed total band energy. Such a resonance implies the existence of two-atom quasi-bound states in the lattice, 
but only if its width (its reciprocal lifetime) is given by the positive value 

Im(D) 
^Re(D) 

Information about the two-atom states can be inferred from Bragg scattering spectra, which depend on the atomic 
density of states in the lattice. The change in the two-atom density of states due to the interatomic interaction 
potential Vab is given by 

^-^M^D)}, (5) 

N being the normalization constant. In the vicinity of a genuine resonance located at E — E r (i.e., if the conditions 
leading to Eq. (4) are satisfied), Eq. (5) reduces to the Lorentzian shape J?J 

T r /2n 

SP = {E-E r f + Yl/A (6) 

This expression measures the enhancement of the two- atom density of states near E r , relative to that of a perfect 
optical lattice (without atom-atom interaction), 

P i0) (E) = -^lm[Tr(G)}. (7) 

Experimentally, however, T r will be swamped by the much larger line broadening due to the Fheat discussed above. 

We now consider that the atoms are located near the minima of the periodic potential wells of the optical lattice, 
for which the tight- binding approximation is valid. As the basis states, we choose the Wannier functions Xn(r — Rp) 
which have the property: J drx n {r — Rp,)HoXn{r — Rp>) — e n {Rp, — Rp,'), £n being the nth band energy eigenvalue 
and Rp the site location. The symmetrized or antisymmetrized product of the Wannier functions of the two atoms is 
denoted by 

|nm; 0) = M[xn{r A - Rp)Xm,{r B - Rp>)± Xn(rB - Rp)Xm(rA - Rp')}, (8) 

where Af is a normalization constant, and + (-) refers to bosons (fermions). Here we assume that /2 ^ //, i.e. no 
lattice-site is occupied by more than one atom. In the plots and estimates given below only bosons are considered. 

In a simple cubic (SC) lattice, Rp — (i/ii + + kjij,)a, where a is the lattice spacing, ^ {i = 1, 2, 3) are integers 
and i,j and k are the unit vectors along the X, Y and Z axes, respectively, we make use of the well-known dispersion 
relation for the tightly-bound energy bands 

€n(k) = A„(0) + 2A„(l)[cos k x a + cos k y a + cos k z a\. (9) 
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Here A n (0) = J d 3 r AX * n {f A - R p )H^ Xn (r A - i?-), A„(l) = J d 3 r AX * n {f A ~ Rp)H^ Xn (r A - R^), and i?^ = R p + ua, 

u being any one of the unit vectors i,j and k. In the harmonic approximation of the lattice potential, Xn (r) is the nth 
level 3D harmonic oscillator eigenfunction. Under these conditions, the unperturbed Green function can be expressed 
in the chosen basis as || 

G nm {p,,il') = (nm;p\ 1 |nm; jf) 

/>oo 

= -iC / dtexp(iE't)[J ai+ (t)J^ 2+ (t)J^ + (t) + J ai _(t)J^ 2 _(t)J U3 _(t)}, E' <3 
Jo 

/•OO 

= C/ dtex P (-^t)[J Ml+ (t)J m (t)J^(t)+7 Ml _(t)J tt ,_(t)J Ms _(t)], £'>3 (10) 



o 

Here \ii± — ± J„ and /„ are the nth order Bessel functions, E' = {E — A„(0) — e m (0)}/A nm (l) are the 
normalized energies and C = 2M 2 / (2ir\ nm (l)), with A nm (l) = 2{A„(1) + e m (l)}. 

Let us take the circularly-polarized additional laser field to have direction cosines (I/a/3, 1/V3, l/v3) in a SC 
optical lattice formed by far-off resonance blue-detuned lasers ||. Then, for two atoms lying along either of the X, Y 
or Z axes, the interatomic potential in Eq. (1) becomes 

8iTa 2 Ik 2 cos(kr/y/3) cos(fcr) 

V A b = 7. ■ (11) 

6c r 

The crucial point is that for two atoms residing at the nearest lattice sites, this potential is attractive for r ~ Tt/k^ <C 
ir/k. By contrast, if the two atoms lie along the direction (0, A=, — ^) , the corresponding potential 



Vab = cos{V2kr/V3) 



5 . fcsin(fcr) cos(fcr) 

■ - cos fcr + ^ + 

3 r z r A 



(12) 



is repulsive for atoms at adjacent sites in this limit. For near-zone separations r ~ Tr/k^ <C 7r/fc, Eq. (11) reduces 
to the "gravity" -like \/r attractive potential. The competition of such a potential with the van der Waals (VdW) 
short-range repulsion (s-wave scattering) is analyzed in M. However, these VdW effects are irrelevant for atoms 
tightly bound at lattice sites separated by 100 nm or more. We therefore can impose a cut-off separation r c (> ro), 
where ro is of the order of 10 3 Angstroms, for calculating the low-energy collisional dynamics of the two atoms. 

We shall illustrate the effect for two bosonic Li atoms that are in the lowest band (n = m = 0) of a SC lattice 
in the absence of the additional laser (Vab = 0). The Lorentzian change in the density of states [Eq. (6)] is 
associated with resonance states formed (for an appropriate Vab) within the energy range of the two unperturbed 
atoms —3 < E' < 3, that is, as long as the two atoms remain confined within the lowest band. In contrast, a genuine 
bound state (a resonance with zero width F,. = in Eq. (4)) is formed at an energy exceeding that of the two 
unperturbed atoms in the lowest band (E' > 3), that is, when at least one atom populates a higher band in the 
presence of Vab ■ 

From Fig. 1 we see that for the laser propagation direction as in Eq. (11), the interatomic potential is attractive, 
only at cos 8 = y/l/3, for r ~ n/k^ -C n/k. Hence, only these atoms would contribute to the signature of a resonance 
(Fig. 2 - inset): the drastic increase of the two-atom density of states, i.e., the appearance of distinct peaks around 
E r (yet bearing in mind that Thoat is their dominant experimental linewidth). The E r and Eb in Fig. 2 are shown 
to scale roughly linearly with laser intensity (at a fixed detuning S). For the 6 Si/2 — > 6 -P3/2 (852.1 nm) transition 
in Cs, the recoil energy En < 10 kHz, whence conditions (a), (b) are best satisfied by a near-UV lattice, ki, ^S> k, 
F e (kr) > 1, /ld(£ ■ r) < 1 and an additional laser with S - 10~ 4 , I < 0.1 W/cm 2 . 

The solution of the Schroedinger equation, H"i> = E^ for the tightly-binding limit and in the presence of Vab can 
be written as 

* = ^2^2c nm \nm;fiv). (13) 

This is a superposition of the two-atom Wannier functions (8) with coefficients C nm obtainable from the Lippman- 
Schwinger matrix equation pj 

^(nrn; fJ,i/\(H - E)\nm; (J,' 'v')C nm (Rfj,, R v ) 

fi'u' 

+ ^ V nm -nv^'u'C mn (R^ 7 R v i) = 0, (14) 
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where Vnm-^v^ 1 v' = (nm; ixv\VAB\nm\ li!v'). Equation(14) can be solved by using the unperturbed Green's function 
Eq.(lO). The solution is given by 

CnmC-R/U) R v ) — C ? ^(i? M , R u ) + G nm (fJ,V, fl'f' 

)^ / nm;/.i / v' ' v" C nrn {Rp, ff , Ry" ) , (1^) 

fi' ' v' ' v" 

where Cnm pertains to the uncoupled superposition of two-atom Wannier functions satisfying the equation 

Y,(nm W v\(H -E)\nm;»,V)C£l(R fi ,,R v ,)=0. (16) 

The formation of a quasibound or resonance state gives rise to quantum translational correlations (spatial entangle- 
ment) of the two atoms. This entanglement is seen to correspond to strong deviations of the C nm in Eq.(15) from the 
uncoupled coefficients Cnm, at E near a resonance. Its experimental manifestation would be the non- separability and 
redistribution of the two-atom density, as compared to the non-entangled (unperturbed) density of atoms at adjacent 
sites (Fig. 3). 

We have concentrated here on the essential effects of laser-coupled atom pairs at adjacent sites of a moderately 
populated lattice, but at higher densities we can obtain multi-atom correlated arrays. A mean-field analysis reveals 
the possibility of a bosonic "supersolid" — a self-organized periodic array of atoms, under the influence of dipole-dipole 
interactions induced by a circularly-polarized plane-wave laser (without an optical lattice) Q. 

To conclude, we have demonstrated the possibility of observing a novel effect: the quasi-binding or binding of cold 
atoms, initially residing at adjacent sites of an optical lattice, by a plane-wave off-resonant laser. Their controllable 
spatial entanglement may be used for matter-wave teleportation |]lfj|| and other quantum information applications. 
Experimentally, such giant quasimolecules can be revealed as satellite lines (with energies depending on the binding 
potential depth) shifted relative to the Bragg scattering spectral lines of a probe laser, analogously to excitonic lines 
in crystals. 

This work has been supported by the German-Israeli Foundation. One of us (BD) thanks S. Giovanazzi for many 
helpful discussions. 
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FIG. 1. The laser-induced off-resonant dipole-dipole potential for two Li atoms in their ground states (2S1/2) as a function 
of kr, for three interat omic axes at different angles 6 with the direction of propagation of the circularly polarized laser: cos 6 = 
(solid line), cos6< = ^/2/3 (dashed line), cos# = (dotted line). The laser frequency is blue-detuned by 300 F from the 

atomic transition 2S 1/2 -> 2P 3/2 (wavelength 670.77 nm and linewidth T/2n = 5.9 MHz), 1 = 5 W/cm 2 , S = 10~ 4 . Inset (top) 
shows the geometry of the system. The other inset shows the x-component of the periodic lattice potential V x = Vbsin 2 fci,:r 
generated by lasers blue-detuned by 10 4 T from the same transition (k ~ Ul) with S = 7.6 x 10~ 5 
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FIG. 2. Main figure: the dimensionless resonance energy E' r = — 

£,/ _ E-2A (0) 



4A (1) 



4Ao^i{°" > ( mar ked by "+") and the bound state energy 
(marked by "*") as a function of laser intensity / for the parameters of Fig. 1. Here Ao(0) = 3.408 MHz, 



Ao(l) = 0.98 Hz. Inset: the drastic increase of the two-atom density of states 8p by 4 orders of magnitude (solid line) compared 
to the unperturbed scaled density of states p = 10 4 x (dashed line) near the resonance E' r = 2, for low laser intensity 
I = 0.48 mW/cm 2 , 8 = 300r. The resonance linewidth T r = 6.29 Hz. 
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FIG. 3. Inset: Scattered intensity |C(/I)| 2 as a function of laser intensity I (S = 300r). The incident atoms (J = 0) 
are assumed to be in sites (000,111). Here solid line refers to |C(000, 111)| 2 , dashed lines to |C(000, 100) | 2 , "+" to 
C(000, 110) | 2 and dotted line to |C(000, 200) | 2 . Main figure: Two-atom wavefunction squared (times the unit cell volume 
a 3 ) |*l/(r)| 2 = § d A R\^>{R,r)\ 2 , where R is their center-of-mass coordinate, as a function of the separation r (in units of a) for 
I = 0.48 mW/cm 2 (dashed), 1 = 5 W/cm 2 (dotted) and 1 = (solid). All adjacent sites become partially occupied by the 
nonseparable (entangled) atoms 
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